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Introduction. There is a classical formula 
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where a(u) and p(n) are the usual elliptic functions. As far as the author knows the earliest 
work in which this formula appeared is a paper of Frobenius and Stickelberger [10]. Before 
publication of their paper, Kiepert [14] gave the formula 
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which also can be obtained by a limiting process from (0.1). These formulae are included 
in the book [21, p. 458 and p. 460] as exercises and are also treated in [9, I, p. 183 and 
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p. 186]. The function o~(nu)/a(u) n in (0.2) has been denoted by ^> n {u) ever since the book 
of Weber [20]. 

There is also a detailed analytic theory of Abelian functions for the Jacobian variety 
of any hyperelliptic curve. In the theory there is a nice generalization by H.F. Baker and 
others of the elliptic sigma function. If the genus of the curve, say C, is g, this function is 
also denoted by cr(u) = a{u\, ■ ■ ■ ,u g ) (see (1.1) below). Although the function 

a(nu) 

a(u) n2 
l 
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is a function on the Jacobian variety and is a natural generalization of i/> n (it) for the given 
genus g curve, it has poles along a theta divisor 6 which is the (g — l)-fold sum of the 
curve C embedded with a Weierstrass point as base point, so does not naturally restrict 
to a function on the curve itself. 

Now we restrict ourselves to consider only functions associated to a hyperelliptic curve 
of genus g = 2. Then we consider the modified function 

a(nu) 
a 2 (u) n2 1 

where 02 (it) = (da / 81x2) (u). Although this function is not a function on the Jacobian 
variety of the curve, via restriction to C it can be regarded as a function on the curve 
itself, and is a good candidate to be called the hyperelliptic psi function. So we denote it 
by V>n(it). The zeroes of ip n (u) on the curve are just the points whose multiplication by n 
are on 0. If g = 2 the theta divisor 6 coincides exactly with C. 

The aim of this paper is to give a natural generalization of the expressions (0.1) and 
(0.2) for any hyperelliptic curve of genus two. See Theorem 2.3 and Theorem 3.3 below. 

The function ip n (u) has been investigated in several contexts as follows. First, it had an 
important role in the new theory of complex multiplication due to D. Grant ([12] and [18]). 
Secondly, a recursion relation for any hyperelliptic psi function was given by D.G. Cantor 
in [7], and such relation is used to compute the torsion points on the corresponding curve. 
He also gave a determinant expression for ip n (u), which is different from our expression 
and should be regarded as a generalization of the formula of Brioschi [6, p. 770, £.3] to all 
hyperelliptic curves. 

The formula (0.1) also has a classical generalization to all non-singular algebraic curves in 
terms of the Klein prime form, which is the formula (44) in [8, p. 33], and a generalization 
to all hyperelliptic curves of genus two in terms of Gunning's prime form, which is the 
formula (2) of [13]. The referee of this paper kindly informed the author about these two 
kinds of generalizations. 

The author would like to acknowledge that the proof of the main theorem (Theorem 
2.3) was quite simplified by the referee. 

The first thing which inspired the author in this work are the very important and strange 
formulae which appeared in the famous papers [15] and [16]. These papers solved the prob- 
lem of the determination of the argument of cubic and quartic Gauss sums, respectively, by 
using the formulae given by Propositions 7.3 and 7.4 in [15, p. 181] and by (6.37) and (6.38) 
in [16], which resemble quite closely (0.2) above. Matthews mentions in [15, p. 179] that 
such a line of attack was suggested by observations of S.J. Patterson. The author hopes 
that our new formulae will be used to establish a new theory of complex multiplication. 

The author was led to believe in the existence of our determinant expressions through 
continual communication with Matsutani. The author found another generalization of the 
formula of Frobenius and Stickelberger to the three variables case (i.e. n = 2) for genus-two 
curves in [5, pp. 96-97] after having already obtained these results. 

Added for the finally revised version. Quite recently all our formulae were generalized to any hyperellip- 
tic curves by the author himself in [19]. Moreover the connection of Theorem 3.3 below and the determinant 
expression of Cantor mentioned above was completely clarified by S. Matsutani (see Appendix of [19]). 
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0. Convention. We denote, as usual, by Z and C the ring of rational integers and 
the field of complex numbers, respectively. In an expression of the Laurent expansion of a 
function, the symbol (d°(zi, z 2 , • • • , z m ) > n) stands for the terms of total degree at least 
n with respect to the given variables Zi, z 2 , • • • , z m . 

In cross references, we indicate a formula as (1.2), and each of Lemmas, Propositions, 
Theorems and Remarks also as 1.2. 



1. The Sigma Function in Genus Two. In this Section we summarize the 
fundamental facts used in Sections 2 and 3. Detailed treatment of these facts are given in 
[1], [2] and [3] (see also Section 1 of [18]). 

Let 

f(x) = A + Xix + X 2 x 2 + A 3 x 3 + A 4 x 4 + A 5 x 5 

be an polynomial of x over C such that its all zeroes are different each other. Let C be a 
smooth projective model of the curve of genus 2 defined by y 2 = f{x). We denote by oo 
the unique point at infinity. In this paper we suppose A5 = 1. The set of forms 

dx xdx 



2y 



2y 



is a basis of the space of differential forms of first kind. Let 

(X 3 x + 2A 4 z 2 + 3X 5 x 3 )dx x 2 dx 



2y '* 2y 

Then r] (1) and r/ (2) are differential forms of the second kind without poles except at 00. We 
fix generators ai, a 2 , Pi, Pi of the fundamental group of C such that their intersections 



are a, • a,- = A • Pj = 0, a; ■ Pj = for i, j 



UJ 
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1, 2. If we set 

1(3, ^ 4 ^2 



the lattice of periods of our Abelian functions appearing below is given by 



A 



UJ 



+ 00 



(c C 2 ). 



The modulus of C is Z :- 



UJ 



-v. 



We also introduce matrices 



f m f 

_ f ai V2 f a 



v 



Let J be the Jacobian variety of the curve C. We identify J with the Picard group 
Pic°(C) of linear equivalence classes of divisors of degree of C. Let Sym 2 (C) be the 
symmetric product of two copies of C. Then we have a birational map 

Sym 2 (C) Pic°(C) = J 

{Pi, P2) ^ the class of Pi + P 2 - 2 • 00. 

We may also identify (the C-rational points of) J with C 2 /A. We denote by k the canonical 
map C 2 — ► C 2 /A and by 1 the embedding of C into J given by mapping P to the class of P— 
co. Then the image i(C) is a theta divisor which is usually denoted by 0. Although we 
treat only the case of genus 2, we use both of these two symbols to make as clear as possible 
the shape of our expected generalization to the case of higher genus. We denote by O the 
origin of J. Obviously A = k _1 (0) = K _1 i(oo). If u € C 2 , we denote by ui and u 2 the 
first and second entries of u. 
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Lemma 1.1. As a subvariety of J, the divisor is non-singular. 
A proof of this fact is seen, for instance, in Lemma 1.7.2 of [18, p. 390]. If we set 

5" = 
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then the sigma function attached to C is defined, as in [2], by 
a(u) = c exp(— -urj' u' lt u) 

■ exp[2vr^T{i *(n + 8")Z{n + 5") + *(n + 5"){u'~ l f u + 5')}}. * 1 ' ' ' 

n€Z 2 

To fix the constant c above we need the following lemma. 

Lemma 1.2. The Taylor expansion of a(u) at u = (0,0) is of the form 

ui + ^2«i 3 - ^A 5 u 2 3 + (d°(u 1 ,u 2 ) > 5) 
up to a multiplicative constant. 

Lemma 1.2 is proved in [4, p. 96] (see also [11, pp. 129-130] or Proposition 2.1.1(2) in [18]). 
We fix the constant c in (1.1) such that the expansion is exactly of the form in 1.2. 

Lemma 1.3. Let t be an element in A. The function u ^ a(u) on C 2 satisfies the 
translational formula 

a(u + £) = x{?)o-(u)ex.pL(u + e,£), 

where %(•£) = ±1 is independent of u, L(u, v) is a form which is bilinear over the field of 
real numbers and C-linear with respect to the first variable u and L{l w , £ {2) ) is 27T\/— T 
times an integer if £ (1) and £ (2) are in A. 

The details of 1.3 are given in [1, p. 286] (see also [18, pp. 395-396]). 

Lemma 1.4. (1) The function a(u) on C 2 vanishes if and only ifue k~ 1 (Q). 

(2) Suppose v w , v {2) are two points of k _1 /,(C). The function u ^ a{u — v (1) — v {2) ) for 
u £ K _1 t(C) is identically zero if and only if v {1) + v (2) is contained in k~ 1 (0)(= A). If 
the function is not identically zero, it vanishes only at u = v m and v (2) modulo A to order 
1 or to order 2 if these two points coincide. 

(3) Let v be a fixed point of k _1 /,(C). There exist a point v {1) of k~ 1 l{C) such that the 
function u i— ► a(u — v — v {1) ) on K _1 t(C) is not identically zero and vanishes at u = v 
modulo A of order 1 . 

The assertions 1.4(1) and (2) are proved in [1], pp. 252-258, for instance. The assertion 
1.4(3) obviously follows from (2). 
We introduce the functions 

p jk (u) = -^T^loga(u), p jk ... r (u) = -^-Pk-Au). 

These functions are periodic with respect to the lattice A by 1.3, and has poles along 
k _1 (B) by 1.4(1). We also use the notation 



The following formula was also obtained by Baker in [2], p. 381 
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Lemma 1.5. We have that 

o~(u + v)o~(u — v) 
a(u) 2 a(v) 2 



Pll(u) - pll(v) + p\2{u)p22{v) ~ Pl2{v)p22{u). 



Let (u\,U2) be an arbitrary point not in k 1 (0). Then we can find a unique pair of 
points (xi,yi) and (3:2,2/2) on C such that 

U\ = / 0J\ + / Wl, U 2 = / UJ 2 + U 2 (1-2) 



00 ./ oo J 00 



with certain choices of the two paths in the integrals. Then one can show that ([2], p. 377) 

P12H = -xix 2 , P22H = x 1 + x 2 . (1.3) 

If u is a point on k~ 1 l(C), the x- and y-coordinates of t _1 fc(ii) will be denoted by x(u) 
and y(ti), respectively. 

Lemma 1.6. If u is a point on k _1 /,((7), then 

Pl2 / \ Ci 

(it) = \U) = —x(u). 

P22 o- 2 

These equalities are shown in [12, p. 124] by using 1.4(1) and (1.3) (see also the proof of 
Proposition 2.1.1 in [18]). By 1.2 and 1.4(1), we see the following. 

Lemma 1.7. If u e k _1 /,(C), then 

ui = \u 2 3 + (d (u 2 ) > 4) 

Hence 112 is a local parameter at (0,0). Lemma 1.8 below gives a local parameter at any 
point of K _1 t(C) except k _1 (0). 

Lemma 1.8. Letv {0> = (vi ^v^) be a fixed point on k~ x l{C) and assume v {0> k _1 (0). 
Let u = (iii,u 2 ) be a variable point on k~ 1 l(C). Then the variable u\ — v[ 0) is a local 
parameter at v (0) along K _1 i(C). In other words, the function u 1— ► u\ — v{ 0) vanishes at 
u = v (0) of order 1. 

Proof. Since 

d(u 2 — v ( n 0) ) du 2 du 2 dx 

— = = = xlu) 

d(ui — v^) du\ dx du\ 

we have 

u 2 - „!•) = x ( v (0) )( Ul - v[ 0) ) + (d°( Ul - v[ 0) ) > 2). 

Lemma 1.4(3) states that there exists a point v (1) in k~ 1 l(C) such that the function u <— > 
a(u — v (0) — v (1) ) is not identically zero and vanishes at u = v (0) of order 1. Then we have 

= a^-v^Xm - <>) + o-2(-v (1) )(u 2 - <>) + (d°(«i - v[ a \u 2 - <>) > 2) 

Hence the vanishing order of u 1— ► u\ — i^ 0) at u = v (0) can not be higher than 1. So we see 
that u\ — v[ 0) is a local parameter at v m . □ 
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Lemma 1.9. (1) Let u be an arbitrary point on Then a 2 (u) is if and only if 

u is a lattice point, that is, the case k{u) = i(oo). 

(2) The Taylor expansion of the function a 2 {u) on K _1 i(C) at u = (0, 0) is of the form 

a 2 ( u ) = -u 2 2 + (d°{u 2 ) > 3). 

(3) Let v {0> = (v { ±\v 2 oy ) be a fixed point on k~ 1 l(C) and assume vo k _1 (0). The Taylor 
expansion of the function u i— > a{u — v (0) ) on k _1 /,((7) at u = (0, 0) is of the form 

a(u - v {0) ) = a 2 (v {0) )u 2 + (d°(u 2 ) > 2). 



Proof If o~ 2 (u) = 0, then the second equality of 1.6 yields o~i(u) = 0. This contradicts 
to 1.4(1), (2) and 1.1. So it must be cr 2 (u) ^ 0. The assertion (2) follows from 1.2 and 1.7. 
Because o~ 2 (u) is an even function, the assertion (3) follows from 1.4(1) and 1.7. □ 

Definition-Proposition 1.10. Let n be a positive integer. If u G K _1 i(C) ; then 

i>n{u) := 2 . 
o- 2 (u) n 

is periodic with respect to A. In other words it is a function on t(C). 

For a proof of this, see [12, p. 124] or Proposition 3.2.2 in [18]. 

Lemma 1.11. We have that ip 2 (u) = 2y{u). 

This is a result stated in [12, p. 128] (see also Lemma 3.2.4 in [18]). We end this Section 
by stating the following easily shown relations (see Lemma 2.3.1 in [18]). 

Lemma 1.12. If u e K _1 i(C) then 

x(u) = \ + (d°(u 2 ) > -1), y(u) = — L + (d> 2 ) > -4). 
U 2 Z u 2 a 



2. A Generalization of the Formula of Frobenius-Stickelberger. Let us start 
by stating the simplest case of our generalization. 

Proposition 2.1. Assume u and v belong to k~ 1 l(C). Then 



a{u + v)a{u — v) 
a 2 (u) 2 a 2 (v) 2 



= —x(u) + x{v) 



1 x(u) 
1 x(v ) 



Proof. We give two proofs. First we use the formula of 1.5. For u ^ k 1 l(C) and 
v $l k _1 z.(C), after dividing the formula of 1.5 by 



P22(u)p 22 (v) 



(o- 2 (u) 2 - a 22 {u)a(u)) (a 2 (v) 2 - a 22 (v)a(v)) 



£7 (it) 2 



a(v)' 



and bringing u and v close to any points on k 1 i(C), we have the desired formula because 
of 1.4(1) and 1.6. 
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Our second proof is done by comparing the zeroes and poles of each side. If we regard 
the two sides as functions of u, they are functions on c(C). We may assume v k _1 (0). 
The zeroes of the two sides coincide and they are at u = v and at u = —v, and both sides 
have their only pole at u = (0,0) of order 2 by 1-9(2) and 1.12. By 1.9(3), the coefficient 
of l/u 2 2 of the left hand side is —1. Thus we have proved the formula once again. □ 

Remark 2.2. Although the formula of 1.5 is a natural generalization of the corre- 
sponding formula for Weierstrass' functions o~(u) and p(u) (Example 1 in [21, p. 451]), the 
formula of 2.1 bears a striking likeness to the formula for elliptic functions. Finding this 
formula was the clue to the discovery of the formula in Theorem 2.3. 

Our generalization of (0.2) is the following formula. 

Theorem 2.3. Let n be a positive integer. Assume u(= u (0) ), u (1) , ■ ■ ■ , u (n) all belong 
to K _1 t(C). Then 

a(u m + u w + ■■■ + u ( ">) Yli<j <r(u w - u U) ) 



a 2 (u^) n+1 a 2 (u^) n+1 ■ ■ ■ a 2 (u^) n+1 



is equal to 

1 x(u (0) ) z 2 (u<°>) y(u<°>) x 3 (u™) yx(u^) ■■■ yx^ n -^/ 2 {u^) x^ n+2 ^ 2 {u^) 
1 x(u (1) ) x 2 (w (1) ) y(u (1) ) x 3 (u w ) yx(u w ) ■■■ yx^-^/ 2 {u w ) x^ n+2 '>/ 2 (u w ) 

1 x(u (,l) ) x 2 (u (n) ) y{u {n) ) x 3 {u {n) ) yx{u {n) ) ■■■ yx^ n ~^/ 2 (u (n) ) x^ n+2 ^ 2 (u < - n '>) 
or 

1 x(u (0) ) x 2 (« (0) ) y(u (0> ) x 3 (u (0) ) yx(u (0) ) ••• x^ n+1 ^ 2 {u^) yxC"- 3 )/ 2 (u (0) ) 

1 x(u (1) ) x 2 {u^) y{u (1) ) x 3 {u w ) yx(u w ) ■■■ x^ +1 ^ 2 {u w ) yx^-^/ 2 {u w ) 

1 x(u {n) ) x 2 (u (n) ) y(u (n) ) x 3 (u (n) ) yx(u (n) ) ••• x (n+1)/2 (u (,l) ) yx {n ~ 3 ">/ 2 {u {n) ) 
according to whether n is even or odd. 

Proof. Proposition 2.1 is the case of n = 1. We prove this formula by induction on n. 
Suppose the points u (0) , u (1) , ■ • • , u (n) do not belong to k~ 1 l(C). We know that each the 
function 

a ( u (0) + „(!) + . . . + „<»)) JT a ( u (i) _ u O)) 
U U) ^ — 

a(uW) n + 1 a(uW) n+1 ■ ■ ■ a{u^) n + 1 

is periodic with respect to A by 1.4(1) and the theorem of square ([17], Coroll.4 in p. 59). 
After multiplying (n + l)-st power of 

P22(n (0) )p 2 2(u (1) ) • • • P22(n ( ">) 



P222(« (0) )P222(W (1) ) • • • P222(^ ( " ) ) 

{g 2 - o 22 )a ^ (on (<r|-<r 2 2)<r {a\ - cr 22 )a ^ (n) ^ 



0-% - 3aa 2 a 22 + a 2 a 222 - 3aa 2 a 22 + a 2 a 222 a% - 3aa 2 a 22 + a 2 a 222 

bringing u U), s close to any points of K _1 i(C), we see by 1.4(1) the left hand side of the 
claimed formula is, as a function of each u u \ a periodic function on K _1 i(C). Now we 
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regard the both sides of the claimed formula as functions on C of the variable u = u {0> , 
and regard the points itj as points on C. We denote the left hand side by fi(u) and the 
right hand side by f2(u). We consider a divisor 

D = (n + 2)oo — u (1) - u (2) u (n) 

of C. Lemmas 1.4(2) and 1.9 show that the divisor (/1) + D is effective, and 1.12 shows 
that the divisor (^2) + D is also effective. Since D is not a canonical divisor and C is of 
genus 2, the Riemann-Roch theorem shows that the space of the functions / such that 
(/) + D is effective is of dimension 1. Therefore the two sides coincide up to a non-zero 
multiplicative constant. We know by 1.4(2), 1.9(2) and 1.9(3) that the coefficient of the 
Laurent expansion at u = (0, 0) with respect to u g of the left hand side is just the left 
hand side of the hypothetical statement of our induction. We also know by 1.12 that the 
coefficient of such the Laurent expansion of the right hand side is just the right hand side 
of the hypothesis. Thus the formula holds. □ 

Remark 2.4. The formula in 2.3 should be regarded as a generalization of the genus-one 
formula 



a(u m + u (1) + • • • + u (n) ) UiKj °(u (t) 



1 x(u {0) ) y(u (oy ) x 2 (u (oy ) yx(u m ) 
1 x(u (1) ) y{u w ) x 2 (u^) yx{u w ) 



x 3 (u m ) 
x 3 (u^) 



1 x(u (n) ) y(u (Tl) ) x 2 (u (n) ) yx(u (nY ) x 3 (u (Tl) ) •• 
rather than of (0.1). Here we suppose the sigma function u{u) is associated with the elliptic 
curve defined by an equation of the form y 2 = x 3 + X2X 2 + X±x + Ao- The size of the matrix 
above is also n + 1 by n + 1. This formula is easily obtained from (0.1) and explains the 
meaning of the multiplicative constant of the left hand side of (0.1). 



3. A Determinant Expression for Generalized Psi-Functions. We give in this 
Section a generalization of the formula of Kiepert [14] (see also Exercise 33 in [21, p. 460] 
or [9, p. 186]). There is a pretty formula: 

Lemma 3.1. Let u and v be belong to k _1 /,(C). Then 

a{u — v) 



lim 

ui— >«i U\ 

Proof. Because of 2.1 we have 

x(u) — x(v) 



1. 



a(u + v) a(u — v) 



vi o- 2 {u) 2 a 2 {v) 2 ui— ui 



Now we bring u\ close to v\. Then the limit of the left hand side is 



lim 

Ml —*V\ 



x(u) — x(v) dx 



(v). 



ui — vi dui 

This is equal to 2y by (1.2). The required formula follows from 1.11. □ 
Remark 3.2. The reader should take care to note that the left hand side in 3.1 is not 

lim ^ (h€K-h(C)). 

2 

Our generalization of the formula (0.2) for \jj n {u){= a(nu)/a2(u) n ) is the following. 
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Theorem 3.3. Let n be an integer greater than 1. Assume that u belongs to k 1 l(C). 
Let j be 1 or 2. Then the following formula holds: 



(1!2!-- 


(n-l)\)^ n (u) 










x'(u) 


(* 2 )'(u) 


y'(u) 


(x 3 )'(u) 


(yx)'(u) 


{x A )'(u) 


x"(u) 




y"(u) 




(yx)"(u) 




x'"{u) 


(* 2 )'» 


y'"(u) 


(x 3 )'"(u) 


(yx)'"(u) 


(* 4 )'» 



(x^-^iu) y^~^{u) {x z y n -^(u) (yxy—^iu) (x 4 )^- 1 ^) 



Here the size of the matrix is n — 1 by n — 1. TTie symbols ' , " , ■ ■ ■ , ( " x) denote 

(\ 2 / \ n— 1 

3=7 J > (dSjJ > respectively. 

Proof. We start with replacing n by n — 1 in 2.3. If we use = x(-u) ^ in the following 
proof for the case of j = 1, we can easily check that the case of j = 2 is also holds. We 
may assume that u ^ (0, 0) modulo A. For a fixed u, we write h = (hi, h 2 ) = u (1) — u and 
varies u (1) . Then hi is a local parameter at the point u by 1.8 and the right hand side of 
2.3 is equal to 




by Talyor's theorem in one variable. After dividing both sides by hi and taking the limit 
when hi — > while keeping u + h € K _1 i(C), by using 3.1, we arrive at the formula 



a(2u + u (2) + • • • + - u (2) ) 2 ■■■a(u- u'- 1 ') 2 n 2 <i<j <^ (i) ~ uU) ) 



a 2 (u) 2n a 2 (u^) n ■ ■ ■ a 2 (u^-^) n 



1 x(u) x 2 (u) 

x'(u) (x 2 )'(u) 

1 x(u (2) ) x 2 (u (2) ) 



Applying Talyor's theorem in one variable again we see that the right hand side with 
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h = u {2> — u is given by 



1 x{u) x 2 (u) 

x'{u) (x 2 )'(u) 

x(u + h) - (x\u) + x'{u)hi) x 2 (u + h) - (x 2 (u) + (x 2 )'(u)/ii) 

1 x(u (3 >) x 2 {u^) 



1 x(u) x 2 (u) 

x'(ii) (x 2 )'{u) 

^x"(u)(fei) 2 + (d°(hi)>3) ^(x 2 )"( U )(/i 1 ) 2 + (d°(/ ll ) >3) 

1 ' x(u (3) ) x 2 (u (3 >) 



After dividing both sides by (hi) 2 and taking the limit when h (1) — ► 0, by applying 3.1 we 
have 



(1!2!) 



<r(3it + ti (3) + • • • + u ( "- 1) )cr(ti - u (3) ) 3 ■ ■ ■ a(u - n'"" 1 )) 3 U3<i<j ~ uU) ) 



a 2 (u) 3n a 2 (u^) n ■ ■ ■ a 2 (u (n -^y 



1 x(u) x 2 (u) 

x'{u) (x 2 )'(u) 

x"{u) {x 2 )"{u) 

1 x(u <3 >) x 2 {u^) 



Using similar operations repeatedly, we have the formula for ip n (u). □ 
Example. We have 

d t \ _ i / du i t \ _ 2 y f \ tf \ _ D f( x ) 



ani I ax x 



x 

2 



*» = Ji- (I) = 2 X - Z3/( x X 3 ) V , (* 2 )» = %(u), (, 2 )» = ±Df{x), 

where f(x) is as in the beginning of Section l and Df(x) = -^f(x). Then —ip^iu) = 
— a(3u)/a 2 (u) 9 = 8y 3 (u) by 3.3. This example is mentioned in the proof of Lemma 2 (e) 
of [12], too, and the proof here is different from the one there (see also Lemma 3.2(2) in 
[18]). 
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